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We investigate the thermodynamics of a pion gas within the O(A^) model in the 1 /N expansion. 
Using the auxiliary field technique, we compute the effective potential up to the next-to-leading 
(~| order (NLO) and show that it can be renormalized in a temperature-independent manner. The 

crucial step for the consistency of the calculation turns out to be the elimination of the auxil- 



iary field prior to renormalization. Subsequently, we solve the NLO gap equation for the chiral 

D 

(— I condensate as a function of temperature both in the chiral limit and with explicit symmetry break- 

ing. We propose a simple semi-analytic estimate of the NLO correction to the condensate and 
^— I compare it to the exact numerical solution. Finally, we show that in the chiral limit the chiral 

symmetry is restored at finite temperature by a second-order phase transition, and determine the 
critical scaling of the order parameter. We study the dependence of the critical temperature on the 
renormalized coupling and find that in contrast to the weak-coupling limit, at strong coupling the 
critical temperature increases at NLO. 
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1. Introduction 



The low-energy physics of Quantum Chromodynamics (QCD) is governed by chiral symme- 
try, associated with the presence of (almost) massless quarks. In the sector of the lightest u and 
d quarks, the spontaneous breaking of the SU(2)l x SU(2)r chiral symmetry in the vacuum gives 
rise to an isospin triplet of pseudo-Nambu-Goldstone bosons, the pions. As the lowest excitations 
of the ordered ground state, these dominate the thermodynamics at low temperatures. Since the 
conventional perturbation theory breaks down at finite temperature, one needs a suitable nonper- 
turbative resummation scheme to deal with the thermodynamics, in particular with the symmetry- 
restoring phase transition. In QCD with two quark flavors, one can make use of the isomorphism 
of the chiral group with the 0(4) rotation group. The pions as well as the sigma meson are de- 
scribed using the scalar 0(4) model, which is subsequently generalized to 0{N), allowing for a 
nonperturbative (one-particle-irreducible) l/N expansion [1]. 

In this contribution we further develop the calculation of [2, 3] at the next-to-leading order 
(NLO) in the 1/A'^ expansion. In particular we show that the NLO effective potential can be renor- 
malized by temperature-independent counterterms. This allows us to renormalize and solve the 
NLO gap equation. From the methodical point of view, our main conclusion is that for a consis- 
tent renormalization one has to dispense with the auxiliary field introduced in [1] (see [4] for an 
alternative approach). From the physical point of view, we show that the NLO truncation predicts a 
second-order chirally-restoring phase transition and the critical behavior in accordance with general 
universality arguments. For details of the calculations as well as a more complete bibliographical 
account of related work we refer the reader to [5]. 



2. The auxiliary field technique 

At zero temperature and density, the 0(A'^) sigma model is defined by the Euclidean Lagrangian 

where i = I, . . . ,N, and the pai^ameters and A denote the pion decay constant and coupling, 
respectively. The //-term accounts for explicit symmetry breaking by nonzero quark mass. The 
essence of the auxiliary-field trick [1] is to add a new field a by means of a pure Gaussian action. 
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In case we wish to study pion gas at finite isospin density and pion condensation [6-8], we further 
have to introduce the isospin chemical potential jXi in terms of the covariant derivative of the scalar 
field. Choosing 2 to represent the charged pion TT^, this amounts to adding 



As a next step one defines the chiral condensate, 0o = {^n)/\N, the pion condensate, po = 
(0i)/\/yV, and the condensate of the auxiliary field, IM^ = {a). The fluctuations of the fields 



2 



0{N) model in 1 /N expansion 



Tomas Brauner 



are suppressed with respect to the condensates by the factor 1 / ^/N. This ensures that the NLO 
contribution to the effective action is given by a simple Gaussian integral over the fluctuations. The 
effective potential up to NLO thus becomes 

1 T T T T NM^ 1 T T 1 _/■ T T 

- IfWiPM)] - ^^In [(p2 + m2)2 + 4M/co2] +^ln(p2+M2), (2.1) 

where = — jij, (On is the bosonic Matsubara frequency, and we used a compact notation for 
the sum- integral over the four-momentum P = {io)„,p). Also, we denoted 

The leading-order (LO) effective potential is given by the first line of Eq. (2.1). The sum-integral 
involved is divergent, but the divergence can, after regularization using a four-dimensional Eu- 
clidean cutoff, be absorbed into renormalization of the parameters /2 and 1/A. The requirement of 
the stationarity of the effective potential leads to a set of gap equations, whose solution yields the 
standai^d phase diagram [6]. 



3. Next-to-leading order 



The NLO effective potential is given by the second line of Eq. (2.1). The last two terms 
represent the correction to the pressure of the free pion gas induced by a finite chemical potential. 
With a four-dimensional cutoff (which is crucial for the successful renormalization at the next-to- 
leading order), they give rise to an unphysical /i/-dependent quadratic divergence. In the following, 
we therefore restrict ourselves to finite temperature, but zero chemical potential. 

In general, the NLO effective potential has to be evaluated numerically. However, its divergent 
part can be calculated analytically [2, 3]. The coefficient of the quadratic divergence depends 
explicitly on temperature; this dependence disappears only when the LO gap equation for M is used. 
This suggests the necessity for a reinterpretation of the effective potential [5]. Since the auxiliary 
field a was introduced as a composite operator of 0,, it may give rise to unphysical features of the 
effective potential. Once it is eliminated using its LO equation of motion, the resulting effective 
potential as a function of 0, alone can be renormalized in the usual manner. Using this elimination 
procedure, we renormalized the NLO effective potential and solved the gap equation for the chiral 
condensate. The results are shown in Fig. 1. The numerical values were obtained with fji = 47 MeV 
(differing by a factor 1 /2 from the conventional value), A = 30, renormalized at the scale 100 MeV, 
and H = (104 MeV)^, adjusted in order to reproduce the physical pion mass in the vacuum. 

Apart from a direct extremization of the effective potential, one may attempt to solve the gap 
equation by a systematic expansion in powers of 1/A'^. Write generally the effective potential per 
degree of freedom as V{^) = Vo{<p) +xVi(0), where Vb.i denote the LO and NLO contributions. 
Appai^ently, setting x = restricts the effective potential to its LO value, while xo = l/4 recovers the 
NLO result. Treating x as a continuous variable, i.e., demanding that the gap equation dV /d<p = 
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Figure 1: Chiral condensate as a function of temperature in the chiral limit (left panel) and at the physical 
point (right panel). LO (dashed line) and NLO (solid red line) values are compared. For explanation of the 
green and blue lines see the text. 

be satisfied for all x G [0,xo], one derives a formal integral equation for the condensate, 



0M = 0(O) 



dx - 



/o v^'{Hx))+xvi'{(^{x)y 

Assuming that the second derivative of the effective potential is constant (that is, approximating it 
at the LO point by a parabola), we find the explicit expression 

This is the green line in Fig. 1. Note that the formal 1/A'^ expansion to the next- to-leading order 
would require neglecting V[' in the denominator (blue line in Fig. 1). Obviously this is not a good 
approximation: It overshoots the NLO coiTcction to the chiral condensate by nearly 100%. 

In Fig. 2 we display the dependence of the NLO critical temperature on the renormalized cou- 
pling. While the weak-coupling analytic calculation predicts a negative correction, Tc = ./ r^fn, 

V 

the numerical results show that at strong couphng, the critical temperature is increased by the NLO 
correction. Even though the auxiliary-field technique does not allow us to check the weak-coupling 
Mmit directly, the decreasing trend towards smaller values of the coupling is reassuring. Finally, 
from the solution of the gap equation near the critical temperature we were able to extract the 
critical exponent, which governs the scaling of the order parameter. The numerical value is in 
agreement with the analytic NLO expression, V = ^ — 

To conclude, we have shown that the NLO effective potential of the 0(A'^) model can be con- 
sistently renormalized in a temperature-independent manner provided one uses the LO equation of 
motion to eliminate the auxiliary field before renormalization. We used this strategy to solve the 
NLO gap equation and proposed a simple semi-analytic approximation to the NLO chiral conden- 
sate which is very close to the exact value. We also investigated the coupling-dependence of the 
NLO critical temperature and found that the NLO correction is negative at weak coupling, but posi- 
tive at strong coupling. In general, the NLO coiTcctions to the observables turn out to be consistent 
with the 1 /N expansion, i.e., roughly of the order of ten percent. 
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Figure 2: NLO critical temperature as a function of the coupling. The dashed line denotes the LO (coupling- 
independent) value, while the dotted line indicates the NLO weak-coupling limit. 
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